Mismatch management for optical and matter-wave quadratic solitons 
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We propose a way to control solitons in x (quadratically-nonlinear) systems by means of periodic 
t modulation imposed on the phase-mismatch parameter ( "mismatch management" , MM) . It may be 

realized in the co-transmission of fundamental-frequency (FF) and second-harmonic (SH) waves in 
a planar optical waveguide via a long-period modulation of the usual quasi-phase-matching pattern 
of ferroelectric domains. The MM may also be implemented by dint of the Feshbach resonance in 
a harmonically-modulated magnetic field in a hybrid atomic-molecular Bose-Einstein condensate 
(BEC), with the atomic and molecular mean fields (MFs) playing the roles of the FF and SH, 
respectively. The problem is analyzed by two methods. First, we identify stability regions for spatial 
^ solitons in the MM system, in terms of the MM amplitude and period, using the MF equations for 

spatially-inhomogeneous configurations. In particular, an instability enclave is found inside the 
stability area.The robustness of the solitons is also tested against variation of the shape of the 
input pulse, and a threshold for the formation of stable solitons is found in terms of its power. 
Interactions between solitons are virtually unaffected by the MM. The second method {parametric 
approximation), going beyond the MF description, is developed for spatially- homogeneous states. 
It demonstrates that the MF description is valid for large modulation periods, while at smaller 
periods the non-MF component acquires gain, which implies destruction of MF under the action of 
the high-frequency MM. 
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CZ) , PACS numbers: 42.65.Tg; 03.75.Lm; 03.75.-b; 05.45.Yv 

I. INTRODUCTION 

Solitons are robust localized pulses that have been predicted theoretically and created experimentally in diverse 
physical settings. The current research in this field is heavily focused on nonlinear optics and Bose-Einstein 
condensation (BEC) 0. Typical solitons are found in uniform media with constant characteristics. However, in 
^vq I many cases it is necessary to consider solitary waves (which are also called "solitons", in a loose sense, even if 
. they are not described by integrable equations) traveling across heterogeneous media, or subjected to strong time 
modulation. A well-known example of the former setting is dispersion management, which is an important concept 
in fiber-optic telecommunications, helping to support stable soliton trains used as data-carrying streams 0,0- On 
the other hand, a possibility to stabilize matter-wave solitons by means of the nonlinearity management, applied to 
CZ2 \ them via the Feshbach resonance in a modulated magnetic field (i.e., time-periodic variation of the scattering length 
that determines the coefficient in front of the cubic term in the Gross-Pitaevskii equation), in one-dimensional (ID) 
and 2D geometries, has drawn considerable attention in the studies of BEC. Another example of a setting 
supporting the transmission of robust solitons in a strongly heterogeneous periodic system, that combines features 
of both the dispersion management and nonlinearity management, is the split-step model (SSM). In the simplest 
case, it is composed of periodically alternating pieces of optical fibers with zero dispersion and zero nonlinearity (i.e., 
it is built as a periodic concatenation of nonlinear and dispersive segments, the latter ones taken with anomalous 
dispersion) 0, 0). In a more realistic variant of the SSM, the nonlinear and dispersive segments are allowed to 
have nonvanishing dispersion and nonlinearity, respectively 0. Multicomponent generalizations of the SSM were 
elaborated too, including one for the WDM (wavelength-division- multiplexed) system 0, and a model taking into 
regard two polarizations of light 0] • 

In the above-mentioned examples, solitons are supported by the cubic nonlinearity of the medium. It is well known 
that the quadratic (second-harmonic-generating, alias x^) nonlinearity also gives rise to stable solitons, that have 
been studied in detail theoretically and experimentally in optics In most cases, these are spatial solitons, i.e., 
self-supporting localized beams in bulk or (which is more relevant to the present work) planar waveguides. Temporal 
v/ 2 ) solitons have been created too 0,0], but under very sophisticated conditions. 

In terms of BEC, a counterpart of the second- harmonic generation is the Feshbach association in atomic BEC 
[TaL IT3 . 0| , induced by the coupling of atomic and molecular mean fields (MFs) with the help of resonant optical 
fields or by hyperfine interactions, using the Zeeman effect for the mismatch tuning. Quadratic solitons in BEC 
have been predicted in Ref. ^(|. In terms of the comparison with optics, these solitons in BEC may be classified as 
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temporal ones. 

In those studies, it was established that conditions for the existence and stability of solitons are most sensitive to 
the mismatch between the fundamental-frequency (FF) wave and the second harmonic (SH), or, in terms of the BEC, 
between the atomic and molecular MFs. This fact suggests a natural question, which is the subject of the present 
paper - how the quadratic solitons will react to a periodic modulation of the mismatch, i.e., "mismatch management" 
(MM). The issue is of general interest, as a possible contribution to the theory of the "soliton management" Q, 
and may also be potentially promising as concerns the use of \^ spatial solitons in optical switching and other 
applications to photonics. As for the soliton-management schemes, almost all of them were explored in terms of media 
with the cubic nonlinearity; the only example dealing with a \^ setting was the model of "tandem solitons" , which 
assumed their transmission in a waveguide built as a concatenation of x^and linear segments, in ID 01 an d 2D [l^ 
geometry (actually, the purpose of the tandem model was to reduce the mismatch). 

A ubiquitous approach to the reduction of mismatch in optical waveguides is based on the use of the quasi-phasc- 
matching (QPM) scheme, in which the material of the x^ waveguide is subjected to periodic poling (since the 
material, such as LiNbC>3, is a ferroelectric, this is usually carried out through periodic reversal of the orientation of 
ferroelectric domains) The poling gives rise to a change of the sign of the x*- 2 ' coefficient with a certain period, 
Lqpm, and thus adds an extra wave vector, Icqpm = (27t/Lqpm) aligned with the propagation direction, e z , to 
the relation between the FF and SH wave vectors, kpp and kgH, which may be used to cancel the original mismatch, 
2kpp — kgn- In terms of this technique, the MM may be implemented by imposing a long-period supermodulation on 
the QPM poling. 

Besides its direct relevance to optics, the MM may also be implemented in the above-mentioned atomic-molecular 
BEC, through the Feshbach-management technique. As mentioned above, in terms of BEC, the application of the 
latter technique, which is based on the Feshbach resonance driven by a modulated magnetic field, to the stabilization 
of various types of ID 2D |(| (and also 3D J2jj) matter- wave solitons was theoretically studied in detail for 
the atomic BEC which obeys the Gross-Pitaevskii equation with the cubic nonlinearity, but no similar results were 
reported, thus far, for atomic-molecular condensates. 

Although in the limit of large values of the mismatch a x system may be reduced to a x*- 3 -* limit by means of 
the well-known cascading approximation there is a fundamental difference between the systems. While x models 
may have exact Bethe-ansatz solutions |21| even beyond the MF approximation, the introduction of the SH (the 
molecular field) lifts the integrability j^. Thus, we expect effects of MM in \^ systems, at small or moderate values 
of the mismatch, to be different from earlier studied effects of the nonlinearity management in the x*' 3 '' model (in the 
cascading limit, the x MM goes over into the x^ 3 ' nonlinearity management). 

The objective of the paper is to study solitons and their stability in one-dimensional MM systems, in both the 
optical and BEC realizations. The MF model, based on a set of partial differential equations, is introduced in Sec. II. 
A piecewise-constant periodic modulation of the mismatch parameter in this model is natural in the optical setting. 
Basic results for solitons in the MF model are collected in Sec. III. We report stability regions for the solitons, and 
conditions necessary for their self-trapping from input beams. Interactions between the solitons are considered too, 
with a conclusion that characteristics of the interactions in the MM system are virtually the same as in its ordinary 
(unmodulated) counterpart. 

In Sec. IV, we employ the parametric approximation |23| . which goes beyond the mean field, and is based on a 
set of ordinary differential equations, in the case of spatially uniform configurations. Although this model is unable 
to g enerate solitons, it allows us to analyze effects of quantum fluctuations, which may be important in both optical 
|24j and, especially, BEC [2^ realizations of the interactions, as well as relaxation effects, which may appear 
in the BEC due to inelastic collisions. In the framework of this analysis, we adopt the usual, in terms of the BEC, 
harmonic form of the periodic modulation of the magnetic field which tunes the Feshbach resonance, rather than the 
piecewise-constant format, adopted in the optical model. Numerical solutions of parametric-approximation equations 
demonstrate that the non-condensate component in the atomic-molecular gas is not essentially excited by the MM, 
under a natural condition that the modulation frequency is low enough. On the other hand, the modulation at higher 
frequencies may lead to destruction of the condensate. The paper is concluded by Sec. V. 



II. THE MISMATCH-MANAGEMENT MODEL FOR THE OPTICAL MEDIUM 

In a normalized form, which is widely adopted in nonlinear optics, the fundamental x*- 2 ^ model in one dimension 
is based on a system of coupled equations for complex local amplitudes of the FF and SH waves, it(z, x) and v(z, x) 
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where a is the mismatch parameter [this coefficient is an irreducible one in the framework of the notation adopted 
in Eqs. QJ, if solitons are intended to be looked for in the z-independent form], and the asterisk stands for the 
complex conjugation. Here, it is assumed that, in the case of spatial solitons, light propagates along axis z in a 
planar waveguide with the transverse coordinate x. In terms of temporal solitons, x is the reduced-time variable. MF 
equations for atomic-molecular BEC in an atomic waveguide, where x is directed along the waveguide axis, can be 
reduced to the form of Eq. JJ) as well, with z being time (see Sec. HVI below^l. 

In the notation of Eq. QJ, complete matching is attained at a = 4, while the single exact analytical solution for 
the x* 2 ' soliton is available at a = 1, 

u = ± (3/V2) sech 2 (x/2), v = (3/2)sech 2 (x/2) (2) 

(the Karamzin-Sukhorukov soliton 26]). At other values of a, solitons were found numerically, as well as in an 
approximate analytical form, by means of the variational method 

The simplest model of the MM (similar, in particular, to that of the dispersion management 3, 4]) assumes periodic 
modulation of a according to the following map (the latter term follows the pattern of the "dispersion-management 
map" Hi): 



, , J a - A«, nL < z < (n + 1/2)1/ , . 

a(z) = < . a / ,/„ u . ' / . , u , n = 0,1,2,3..., (3) 
v ' I ao + Aa, (n + l/2)L<z < (n + 1)L y ' 

where ao is the average value of a, while Aa and L are the amplitude and period of the periodic management. We 
here assume equal lengths, L/2, of the two segments forming the MM cell. Simulations supplementing those reported 
below demonstrate that a change of the relative length of the two segments (with the respective change of the local 
values of the mismatch in them) produces little effect on eventual results, quite similar to what is known about the 
dispersion management 

While the MM map in the piecewise-constant form of Eq. J2J) is most natural in terms of optical waveguides, for 
the application to the atomic-molecular BEC a more natural choice is, as mentioned above, the harmonic modulation 
of the mismatch, corresponding to the periodic time dependence of the magnetic field tuning the Feshbach resonance, 
see Eq. (|14f) below. The experience gained in the studies of various models of the dispersion management Q suggests 
that the piecewise-constant and harmonic formats of the modulation cannot lead to qualitatively different results. 

Equations Q conserve a known integral of motion, viz., Manley-Rowe invariant, 



'MR 



\u(x)\ 2 + 4\v(x)\ 2 dx, (4) 



even in the case of the z-dependent mismatch parameter, whereas the Hamiltonian corresponding to Eqs. Q is not 
conserved. In addition, in the case of "walking" solitons [2£j [in fact, these are spatial solitons tilted in the plane of 
(x, z)], Eqs. with variable a(z) conserve the total momentum, 

/+00 
«u + 2v*v) dx. (5) 
-00 

It is also worthy to note that Eqs. QJ with a = a(z) are invariant with respect to the spatial Galilean boost, hence a 
generic tilted soliton, {u, v) , can be generated from the straight one by means of the corresponding transformation, 

u{x, z) = e< c2/i ) z+l(c/2)x u(x - cz), v{x, z) = e< c2 /2 ) z+lcx v(x - cz), 

with an arbitrary real tilt parameter c. As follows from Eqs. iQl - ©, the momentum of the tilted soliton is P = cl, 
i.e., c and I play the role of the effective velocity and mass of the soliton. 

In the case of the atomic-molecular BEC mixture, the Galilean invariancc has its literal meaning (in the temporal, 
rather than spatial, domain). However, it is broken if the axial trapping potential is taken into regard. 



III. DYNAMICS OF SPATIAL SOLITONS UNDER THE MISMATCH MANAGEMENT IN THE 

OPTICAL MODEL 

A. Formation of stable mismatch-managed solitons 

To simulate the transmission of spatial optical solitons under the MM conditions, we solved equations Q with MM 
map @ by means of the split-step numerical method, which uses the Fourier transform to handle the linear stage of 
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(a) (b) 

FIG. 1: (Color online) A typical example of the quick self-trapping of an initial beam into a stable spatial soliton in the 
mismatch-management model, with ao = L = Aa = 1. The input was the Karamzin-Sukhorukov soliton |2(l for a(z) = 1, 
taken as per Eq. Panels (a) and (b) show the evolution of the FF and SH fields. 

the evolution. As the input (initial pulse), we took either an ordinary soliton corresponding to the averaged version of 
the model, i.e., one with a(z) = ao, or a deliberately altered pulse, to verify whether the MM system will provide for 
self-trapping into a stable transmission regime. Below, we display systematic results obtained for ao = 1 and ao = 2 
(comparison with results collected with other values of the average mismatch parameter demonstrate that these two 
cases adequately represent the generic situation). In particular, for ao = 1 we launched the initial pulse taken as 
exact solution (J2J) corresponding to a — 1. For ao = 2, we typically used either a soliton solution found in a numerical 
form for a = 2, or, in order to try the effect of a strong change of the input, we again took expression (J2J), i.e., the 
exact soliton appertaining to a = 1. At a fixed value of ao, results were collected by varying the MM amplitude Aa 
and period L. 

First, in Fig. ^ we display a typical numerical solution for ao = 1, L = 1 and Aa = 1, generated with the use of 
exact soliton J5J as the initial condition. In this figure, panels (a) and (b) present the evolution of the FF and SH 
components of the field. In fact, in all cases considered, there was no conspicuous difference in the dynamics of the two 
components, therefore in other cases shown below we only display the picture for the FF beam. As is seen from Fig. 
n the input beam readily gives rise to a robust spatial soliton (intrinsic pulsations with period L = 1, caused by the 
MM, are almost invisible in Fig. QJ; the soliton remains stable in indefinitely long simulations. The transient stage, 
necessary for the self-trapping, is fairly short, comprising a few MM cells. In the subsequent evolution, gradually 
fading residual oscillations of the pulse's amplitude can be seen, with a period covering several cells (these oscillations 
are caused by the initial perturbation, rather than the periodic MM). Further examples of stable transmission regimes 
are displayed below in Figs. [Sfa) and (c). 

B. Stability diagrams 

Conclusions drawn from systematic simulations are summarized in the stability diagrams, which are displayed in 
Fig. [21 for ao = 1 (a) and ao — 2 (b). As said above, in the former case we launched the pulse corresponding to 
Karamzin-Sukhorukov soliton J2Jl, which is an exact solution for a(z) = 1, and in the latter case the initial pulse 
was a numerically found stationary soliton corresponding to a(z) = 2. The diagrams display areas in parameter 
plane (Aa, L) where the initial pulse gives rise to stable transmission, or decay of the pulse [an example of the latter 
outcome is displayed below in Fig. ffib)] . Naturally, the stability regions tend to extend along the parameter axes, as, 
in either limit of Aa — > or L — > 0, the model returns to the usual x^ system (in the case of L — > 0, this is provided 
by averaging), where the initial pulse represents an ordinary stable soliton. 

Generally, the instability of the x^ solitons in a part of the parameter space may be realized as a result of a 
resonance between the perturbation frequency, introduced by the periodic action of the MM, and the frequency of the 
intrinsic mode, which, as is well known, x^ solitons have in the system with constant coefficients [on the contrary to 
the nonlinear-Schrodinger (NLS) solitons] Indeed, comparison of the numerically found instability border, shown 
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(a) (b) 

FIG. 2: Stability and instability areas in the mismatch-management model, in the plane of the management amplitude (Aa) 
and period (L). Panels (a) and (b) display the stability diagrams for the model with average mismatch ao = 1 and ao = 2, 
respectively. In either case, the diagram was built by collecting results of many simulations, with the input taken as a soliton 
of the respective averaged model, i.e., one with a(z) = ao [for ao = 1, the initial pulse is the analytical Karamzin-Sukhorukov 
solution, given by Eq. J5J]. 



in Figs. Ufa) and (b), with values of the intrinsic eigenfrequency, which are known from numerical computations too, 
demonstrates that the explanation of the instability by the resonance is feasible, although the intrinsic frequency is 
defined for infinitely small perturbations, while the instability sets in when the solitons are strongly perturbed. 

It is also relevant to note that, in the limit of large positive values of mismatch a, the SH (or the molecular MF, 
in terms of BEC) can be eliminated from Eqs. JQ) by means of the well-known cascading approximation, v ss u 2 /2a, 
reducing the system to the single NLS equation. The same is valid when the large mismatch is subjected to the 
management [provided that a(z) does not change its sign]. In that case, the cascading limit will lead (as mentioned 
above) to the NLS equation featuring periodic nonlinearity management 4]. In fact, solitons in the latter equation 
and their stability were investigated in some detail in various contexts 0, U- ^ n particular, a resonant mechanism of 
the destabilization of higher-order solitons (bound states of fundamental solitons) under the action of the nonlinearity 
management, qualitatively similar to one outlined above, has been demonstrated in Ref. |2?| . 

A noteworthy feature of the diagrams is the presence of an instability enclave inside the stability area. In Fig. Ufa), 
the enclave is shown symbolically by a square, as exact delineation of its borders requires extremely long simulations. 
In Fig. |2Ib), the borders of the enclave approximately correspond to its real shape. It may be relevant to mention 
that examination of similar stability diagrams in the above-mentioned SSM had revealed a different but somewhat 
similar feature, viz., a system of stability islands inside the instability area at large values of L Q. In the present 
MM model, the simulations do not reveal stability islands (on the other hand, no "instability lakes" , that would be 
similar to the enclaves in Fig. |2 were found inside the stability area in the SSM). We surmise that there may exist 
additional small instability enclaves, and, in principle, they may even form a fractal pattern. However, an accurate 
investigation of these issues requires extremely high numerical accuracy, and they are left beyond the scope of this 
work. 

The existence of the instability enclave is further illustrated by a set of simulations presented in Fig. |3 which are 
performed along a vertical line, Aa = 0.3, cutting through the enclave in Fig. |2a). A notable difference is observed in 
the self-trapping into stable transmission regimes below and above the instability region: in the former case, at L = 4 
(which is close to the instability border), the established pulse is very different from the input, due to considerable 
radiation loss in the process of the establishment of the stable-transmission regime (in the SH component of the wave 
field, which is not shown in Fig. |3 an approximately the same degree of the loss is observed). It is relevant to mention 
that the formation of stable solitons in the SSM may also be accompanied by strong losses, depending on parameters 
of the system and the form of the input 7]. On the other hand, above the instability enclave, the loss is small, and 
the established soliton is closer to the input, as seen in Fig. [3Jc) (the same is observed in the SH counterpart of the 
latter figure, which is not shown here). 
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(b) (c) 

FIG. 3: (Color online) Three examples of the evolution generated by the input in the form of Eq. J2J launched into the 
mismatch- management system with ao = 1 and Aa = 0.3, the modulation period being L = 4 (a), 6 (b), and 7 (c). Point (b) 
falls into the instability enclave in Fig. |2Ia); in this case, the soliton suffers, eventually, complete destruction. 

C. The system's tolerance to variations of the input pulse 

Another noteworthy feature revealed by the simulations is great tolerance of the spatial solitons established in the 
MM system to variance of the input beam. In particular, if the input launched into the system with ao = 2 was 
deliberately taken in the "wrong" form of expressions (0) (recall they would yield an exact solution for the averaged 
equations with a = 1, rather than a = 2), the simulations, performed for various values of Aa and L, produced a 
stability diagram nearly identical to that shown in Fig. Gib). 

It is natural to expect that decrease of the input power will eventually lead to a failure in the self-trapping of the 
spatial soliton, i.e., there must exist a certain power threshold for the soliton formation. To find it, we performed 
additional simulations with the inputs similar to those used above [in particular, waveform @ was taken for ao = 1], 
but multiplied by a power-reducing factor, W < 1: 

{uq(x), v (x)} -> Vw {uo(a;), v Q (x)} . (6) 

It was found that the critical value of W, below which the thus altered input pulse fails to generate a stable soliton 
is Wcr ~ 0.72 (in fact, it is close to a critical value that can be found numerically in the ordinary system, with 
Aa = 0). Figures E|a,b) show what happens when W is taken, respectively, above and below W CI . 
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(a) (b) 

FIG. 4: (Color online) (a) Self-trapping into a stable mismatch-managed spatial soliton of the input taken as per Eq. @ with 
W — 0.81. (b) Decay of the input beam taken as in Eq. ©, but with W = 0.64. In either case, uo(x) and Vo(x) are components 
of soliton solution @, and parameters are ceo = Aa = L = 1. 




FIG. 5: Interaction of two identical in-phase solitons in the ordinary x^ model (with a = 1, Aa = 0) (a), and its mismatch- 
managed counterpart with ao = Aa = L = 1 (b). The initial separation between the solitons is Aa; = 3.8, while the minimum 
separation at which the solitons do not interact is Aa; m in « 5, in this case. 

It was also checked that the change of the distribution of the Manley-Rowe invariant [see Eq. between the 
two components of the input beam virtually does not affect the self-trapping threshold (in particular, one may start 
with the entire power put in the FF field, while vq — 0). The ordinary x^ system, with Aa = 0, features a similar 
property [TT| . 

D. Soliton-soliton interactions 

We have also performed systematic simulations of interactions between solitons in the MM model. It was found that 
the character of the interaction remains virtually the same as in the ordinary x model with constant mismatch, that 
was studied in detail in earlier works The main characteristic of the interaction is a minimum initial separation 
between co-propagating identical in-phase solitons, Aa; m i n , which is defined so that the solitons do not demonstrate 
any interaction for Aa; > Aa; m i n , while, being placed at distance Aa; < Aa; m i n , they start to attract each other and 
eventually merge into a single beam. A typical example displayed in Fig. [S] shows that the interaction indeed seems 
identical in the MM system and its ordinary counterpart. 
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IV. EFFECTS OF RELAXATION AND QUANTUM FLUCTUATIONS 



As explained in Introduction, the \^ model can also describe a hybrid atom-molecule BEC Il3| under quasi-lD 
confinement. In the parametric approximation, developed for the case of free space in Ref. |23j . the condensate 
is described by the following system of equations of motion for the molecular mean field ip m (X, t) and the atomic 
annihilation operators ^ a (X, i): 



d ,,, , 

l-J-ifm {X,t) 

at 



l J* a (X,t) 



1 d 2 

^dX 2 



Vm {X,t)+g(V a (X,t)V a (X,t)) 



^ ($+ ( X , t) * Q (X, t)) + k m \y m (X, t) | 2 ) Vm (X, t) . 



id 2 l . , 

2moX i 2 



+2g*<p m {X, t) *t (A, t) - -ik a \<p m (X, t) | 2 * a {X, t) + iF (X, t) , 



(7) 



(8) 



where m is the atomic mass and units with fi = 1 are used. As shown in Ref. J22J, the ID coupling constant, g, and 
(time-dependent) detuning Did can be expressed, as follows, in terms of elastic scattering length a^ g , phenomenolog- 
ical Feshbach-resonance strength A, the difference between the magnetic momenta of an atomic pair in the open and 
closed channels, (i, detuning of the external variable magnetic field, B(t), from its resonant value Bo, and transverse 
trap frequency u>x : 



>b g M|A, 



D 1D (t) =/i[B (t) - Bo 



■ U)± 



(9) 



The effect of the confinement-induced resonance [2^| is neglected in Eq. |JSJ, as db g is much smaller than the transverse 
size of the trap. The source of the quantum noise in Eq. JSJ, F(X,t), and rate coefficient k a account for the 
deactivation in atom-molecule inelastic collisions |23j |. Rate coefficients k a and k m (for molecule-molecule collisions) 
are related to their 3D counterparts, k a , m — (mLu±/2ir) ka,m ■ 

Neglecting the collision-induced deactivating, and replacing the atomic-field operator ty a (X, t) by c-number mean 
field tpo (X,t), the substitution 



tpa(X,t) = $ucxp^-i J Q(t)dtj, <p m (X,t) = -V2$vexp f-2i J Q 



(t) dt , 



t = z 



I (2\/2> 



(10) 



X = xj (2 5/4 v /mg$) , D XD = \f2g§ (a - 4) 



with SI (t) = —Did (t) /2 — 2v / 2.g < I > , casts Eqs. Q and JHJ precisely in the form of Eqs. JU, which are adopted in 
nonlinear optics. Necessary normalization of u and v can be provided by the proper choice of the mean-field scaling 
constant, <&, in Eqs. (JTDJ. 

Effects of deactivation and quantum fluctuations neglected in normalized equations Q can be taken into regard 
in the framework of the parametric approximation |23T | , where the A-dependence of ip. m is neglected, and the atomic 
field operator is represented as 

(11) 



9 a (X,t) = (l/v^F) / dpe ipx C(t) [A(p,t)^ c {p,t)+ti{-p,t)^ s {p,t) 
C(t) = exp |-i J dt'k a \ip m {t') 



o 

Here, c-number functions V'c.s (j>, t) satisfy the time-evolution equations 

- „2 ! 



iipc,s (p, t) 



£r-§^ (*) 

2m 2 



ipc,a (p, t) + 2g*ip m (t) c (p, t) . 



(12) 



In the present analysis, the initial moment, t = 0, corresponds to a relatively small detuning. As the initial state 
is implied to be a stable condensate, Eq. (|llfl may be considered, at t = 0, as the Bogoliubov transformation, with 
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operators A (p, 0) being annihilation operators of the Bogoliubov quasiparticles. In this case, initial conditions for 
functions ip c ,s(t) are produced by the Bogoliubov transform, 

1/2 

4>c (p, 0) = (^E+lE) , ^ (p, o) = _ 2 ^W ^ (p, o) 



(for p 7^ 0), where d p = p 2 / (2m) — Die (0) /2 — E 1 , the Bogoliubov excitation energy is e p = d 2 — 4\gip m (0) | 2 , 
and _E is the chemical potential of the atom-molecule condensate. At zero temperature, one has A(p, 0) |in) = 

1/2 

(27r) v Vo (0)5 (p) |in), where | in) is the initial-state vector, and the atomic-condensate mean field <po (t) is expressed 
in terms of solutions to Eq. 1|12[| . which satisfy initial conditions ?/> c (0, 0) = 1, ip s (0, 0) = 0, as follows: 

ipo (t) = (in|* „ (X, t) |in> - C (t) [0 C (0, t) (0) + V, (0, t) <p* (0)] . 
Then, the normal and anomalous densities in Eq. Q become X-independent, as shown in Ref. [23^ : 

f)*«(X,t)) = |v?o (t)\ 2 + ^- I dpn s (p,t) 



2tt . 

— OO 

(i>l(X,t)y a (X,t)) = Vl(t) + ^ J dpn s (p,t), 

— oo 

where the momentum distributions of non-condensate atoms, n s (p,t) and m s (p,t), can be expressed in terms of 
-0c, s (p, t) Accordingly, the equation for the molecular mean field, Eq. Q, takes the form of 

i<Pm (t) = g<fl (t) + - I dpm s (p, t) 

7T 

2 

-i \<po (t) | 2 + ^ J dpn s (p,t) + k m \<p m (i) | 2 I <p m (t) . (13) 

Unlike the 3D case j23|, the ID problem does not require renormalization, as the integral of m s is free of the 
ultraviolet divergence. Nevertheless, it now diverges at zero momentum. This formal infrared divergence is related to 
phase fluctuations and to the absence of true condensate in the infinite ID system. However, in the present work we 
actually consider coordinate-dependent solitons, while the divergence is a consequence of the neglect of the coordinate 
dependence in the parametric approximation. A more careful analysis of the inhomogencous case yields an asymptotic 
estimate, m s (p, i) ~ p 2 for p <C Pmim where the characteristic momentum is inversely proportional to the soliton's 
size, Pmin = l/\/mgQ. Thus, p Illm may be naturally chosen as the lower-integration limit in Eq. Ijl3|l . The resultant 
value of the integral features a weak dependence on the lower limit, as the divergence is logarithmic. 

The characteristic kinetic energy of the atoms, both in the condensate and not belonging to it, can be expressed in 
terms of the characteristic momentum as p^J (2m). The system may be considered as effectively one-dimensional 
if this energy is much smaller than the transverse excitation energy, u>±, i.e., oj\ 3> 27mab g ^A/m, where the total 
initial density of atoms, n = muj± (\(po (0) | 2 + 2\(p m (0) | 2 ) / (2ir), is proportional to Manley-Rowe invariant |@J, in 
terms of Eqs. QJ. 

Further calculations involve a numerical solution of Eqs. (|12l) on a grid of values of p, combined with Eq. Ijl3|l . 
Figure presents the results for the 853 G Feshbach resonance in the condensate of 23 Na atoms, with A = 0.01 G, 
fl bg = 3.4 nm, and [i — 3.65/is (see Ref. 23]), with n — 2 x 10 11 cm 3 /s and lj±_ = 1 x 2n KHz. The deactivation rate 
coefficients, k a — 5.5 x 10 _11 cm 3 /s and k m = 5.1 x 10 -11 cm 3 /s, were taken from Ref. 30]. The detuning D\r> is 
defined as per Eq. Q1U|). with the harmonically modulated mismatch parameter, 

a = a + Aa cos (2nz/L) , (14) 

cf. the modulation map in the optical model given by Eq. ©. The choice of the scaling factor in Eqs. (|10fl which 
corresponds to the Karamzin-Sukhorukov soliton, see Eqs. 101, is $ = y/8nn/ (27mui±). The value of cxq = 0.992 
is chosen so as to make the initial molecular fraction equal to half the total population, as at the center of the 
Karamzin-Sukhorukov soliton. 
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FIG. 6: (Color online) Densities of the molecular and atomic fractions in the condensate (solid and dashed lines), and of 
non-condensate atoms (dot-dashed lines), calculated with Aa = 1 for L = 3 (a) and L = 2 (b), see Eq. 1141 . 

FigureEJa) demonstrates that, for large modulation periods in Eq. (|14f) . the non-condensate fraction remains below 
the level of 10%, and the lifetime due to the deactivation is large enough to allow experimental observation of the 
soliton dynamics described in the previous sections, in terms of the optical medium, in the hybrid atom- molecule 
BEC too. 

However, Fig. b) demonstrates that the non-condensate fraction of the atomic population acquires substantial 
gain for smaller modulation periods. Therefore, the mean-field approach may not be applicable in this region for the 
description of the atom-molecular quantum gas, and new instability mechanisms can be expected. 

V. CONCLUSIONS 

In this work, we have proposed a model of the second-harmonic-generating (x^ 2 -*) system with the mismatch 
parameter subjected to the periodic modulation ("mismatch management", MM). The system may be implemented 
in two altogether different physical contexts: the co-propagation of the FF (fundamental-frequency) and SH (second- 
harmonic) waves in a planar optical waveguide, and atomic-molecular mixtures in the BEC (in the latter setting, the 
atomic and molecular mean fields play the roles of the FF and SH components, respectively). The most physically 
relevant approach to the realization of the MM in these media are, respectively, a long-period supermodulation imposed 
on top of the quasi-phase-matching periodic arrangement of ferroelectric domains in the x*- 2 ) optical waveguide, such as 
LiNbC>3, and the Feshbach resonance tuned by a modulated magnetic field in the atomic-molecular BEC. Accordingly, 
the natural form of the periodic modulation is piecewise-constant in the former case, and harmonic in the latter one. 

The main issue considered in the framework of the mean-field approach was identification of the stability region for 
the mismatch-managed spatial solitons in the plane of two control parameters, Aa and L (the MM amplitude and 
period). In particular, a notable feature of the stability area is the existence of an instability enclave embedded in it. 
Also investigated was the robustness of the solitons against variation of the shape of the input beam, and reduction 
of its power. It was found that the stability of the established regime virtually does not depend on the particular 
shape of the input, as well as on distribution of the total power (Manley-Rowe invariant) between the FF and SH 
component in the input. On the other hand, reduction of the initial peak power by a factor of W < 1 reveals the 
existence of a threshold, W cr ~ 0.72, below which the initial pulse decays. 

In the model of the atomic- molecular BEC, we have demonstrated that the time-evolution equations for the mean 
fields are tantamount to the spatial-evolution equations in the optical model. On the other hand, important issues in 
the context of BEC are the stability of the condensate against generation of fluctuational (non-condensate) components 
in the degenerate quantum gas, and losses due to inelastic collisions. These effects were analyzed within the parametric 
approximation, which goes beyond the mean-field description. Numerical calculations have demonstrated that, quite 
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naturally, the condensate is effectively stable against the periodic perturbations introduced by the MM in the low- 
frequency modulation format, and unstable in the high-frequency regime. 
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